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Abstract 

By using the Redfield form of the master equation, we investigate the decoherence times of a 
double quantum dot charge qubit (DQDCQ) in three different cases, namely when it is coupled 
to (I) the piezoelectric coupling phonon bath (PCPB), (II) the deformation coupling phonon bath 
(DCPB), and (III) the Ohmic bath. It is found that our results for case (I) and (II) are in the 
same magnitude with those obtained via the exact path integral methods, while for case (III), the 
decoherence time is in well agreement with the experimental value. 
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I. INTRODUCTION 

The double quantum dot (DQD) charge qubit [1-5] is one of the qubits that are considered 
to be promising candidates for the reahzation of the building blocks of quantum information 
processing. Its two low-energy states are denoted as the local states |0) and |1), which could 
be controlled via external voltage sources. There exist some effective schemes to prepare the 
initial states and readout the final states of the qubit [6] . It is also known that this kind of 
qubit is coupled to its environment inevitably. Therefore, it is believed that the decoherence 
might be the central impediment for the qubit to be taken as the cell of quantum computer. 
Hence, finding out the primary origin or dominating mechanism of decoherence for the qubit 
is a basic task to overcome the difficulty. 

The recently experimental realization of the coherent manipulation of electronic states 
in a double-dot system [7, 8], which was implemented in a GaAs/AlGaAs heterostructure 
containing a two-dimensional electron gas, stimulated a lot of theoretical interests. Various 
methods have been tried to study the system. For instance, through density matrix simula- 
tion, Fujisawa et al. [9] tried to explain its transportation property. The Born-Markov-type 
electron-phonon decoherence at large times due to spontaneous phonon emission of the 
quantum dot charge qubits was investigated by Fedichkin et al. [10]. In 2005, Vorojtsov et 
al. [11] studied the decoherence of the DQD charge qubit by employing the Born-Markov 
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approximation. Based on a unitary transformation, Wu et al. [12] investigated the decoher- 
ence in terms of a perturbative treatment. Thorwart et al. [13] investigated the decoherence 
of the DQD charge qubit in a longer time with a numerical exact iterative quasiadiabatic 
propagator path integral (QUAPI) method [14], while Liang [15] used an iterative tensor 
multiplication (ITM) method [14] derived from the QUAPI to study decoherence of a DQD 
charge qubit in both the piezoelectric coupling phonon bath (PCPB) and the deformation 
coupling phonon bath (DCPB). They found that the decoherence times of the qubit are 
shorter than the reported experimental ones when the qubit in PCPB and DCPB. 

In this paper, we investigate the decoherence times of the (DQD) charge qubit coupled 
to PCPB, DCPB and the Ohmic bath with another method, through solving the master 
equations. The model Hamiltonian and the spectral functions for the different baths are 
introduced in section II. We then develop the Redfield form of the master equation and 
obtain the decoherence times in section III. Conclusions are given in the last section. 

II. QDQ CHARGE QUBIT MODEL 

The DQD charge qubit consists of left and right dots connected through an interdot 
tunneling barrier. Due to the Coulomb blockade, at most one excess electron is allowed 
to occupy the left and right dot, which defines two basis vectors |0) and |1). The energy 
difference between these two states can be controlled by the source-drain voltage. Neglecting 
the higher order tunneling between leads and the dots, the effective Hamiltonian in the 
manipulation process reads [12, 13] 

^eff = Hs + Hb + HsB, 

where Hs is the Hamiltonian of the QDQ charge qubit, Hb is the Hamiltonian for the 
phonon bath, and Hsb describes the electron-phonon interaction. More explicitly, we have 

Hs = hTcO-^, 

Hb = h^ujqh\hq, 
q 

HsB = ha,J2^Mqbl + M*bg). (1) 



Here, T^ is the interdot tunneling, a^ and (T^ are the Pauh matrix, 6j (bg) are the creation (an- 



nihilation) operators of phonons, fkUg is the energy of the phonons, and Mg = Cg/y^2mqUJqh, 
where Cg are the classical coupling constants of the qubit-phonons system. Having written 
above, we now introduce the spectral density which fully describes the effects of the phonon 
bath [16, 17] 

J{uj)=J2\M/S{u-ujg). (2) 



g 
According to Ref. [12], the spectral density of PCPB is 

■ UJd . f uj\ 
1 sm — 
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^-uj^/2uj} 



(3) 



where ujd = s/d and uoi = s/l, d denotes the center-to-center distance of two dots, / the dot 
size, s the sound velocity in the crystal, and 

_ M f6_ 1 _8_\ 
^'''~ n^gs^ [35^ x35j- ^ ^ 

As usual, M is the piezoconstant, g is the density of the crystal, and x is the rate of transverse 
to the longitudinal of sound velocity in the crystal (see for example Refs. [12] and [13]). As 
in the GaAs crystal s ^ 5 x 10^ (m/s). With the parameters of GaAs [18], Wu et al. [12] 
estimated that Qpz ~ 0.035 (ps)^ . The spectral density of DCPB is also obtained. It is 



J''f{uj)=gdfUJ^ 
with 



UJd . [ UJ 

— sm — 

to \UJd 



■^2 
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(5) 



where S is the deformation potential. In the same paper, Wu et al. also propose a value 
Qdf ~ 0.029 (ps)~^. With the help of the definite spectral density functions of the baths, one 
can investigate the dynamics and then the decoherence of the open qubit. 

We also consider explicitly the case in which the behavior of the original spectral function 
J{ijj) has a simple power-law form for u < uJc- 

J{w) = r]uj'e~'^^'^% T] = const, (6) 

with dimensionless damping strength 77 and a cutoff frequency Uc- On general grounds, the 
linear low frequency behavior of J{w) is expected in basically all condensed-phase electron 
transfer (ET) reactions [16], and the frequency Uc then corresponds to some dominant bath 
mode. For its analytic advantages, we use the Ohmic spectral density [17] by setting s = 1 
in Eq. (6). We also study the case when r] = 0.04, Uc = 0.05 (ps)^"*^ [19]. 



III. DECOHERENCE TIMES 

In this section, we first review the Redfield form of the master equation. In an open 
system, the relevant system is characterized by the reduced density-matrix (RDM) p which 
is defined as a trace over all bath variables of the full density matrix pj. 

P = tTB {Pt) ■ (7) 

The equation for the time evolution of p can be consequently obtained from the Liouville 
equation for prp 

^='j,trs[pAt),H]. (8) 

The Redfield form of the master equation for the RDM is obtained from Eq. (8) by imple- 
menting a series of perturbative approximations [20, 21], namely, the system-bath coupling 
is treated perturbatively up to the second order. The bath is assumed to remain in equilib- 
rium, and then Markov approximation gets involved. In the eigenstate representation of the 

system 

H, \p) = Es l/i) , (9) 

the Redfield equation for the RDM reads [22] 

dpaui^) 



iiv 

— —tiO 



dt 



^^m^Pm^(^) + Yl ^m^«aP«a(^), (10) 

where u^^, = {E^ — Ei,)/h and Rf^u^x is the relaxation or Redfield tensor. The first term on 
the right-hand side of Eq. (10) describes the isolated system evolution, while the second one 
represents its interaction with the dissipative environment. The Redfield tensor describing 
the system relaxation can be expressed as 

a a 

with 

1 f°° 
rt,. = ^y^ dt{{\\HsB{t)\v){p\HsB\n))Be-'--^\ (12) 

^~xu,. = j-,j^ dt{{X\HsB\v){p\HsB{t)\^))Be-'-^-\ (13) 

Hssit) = e'^«*/^iJ5Be-'^^*/^ (14) 

where {■■■)b denotes the thermal average over the bath. 



For the Hamiltonian defined above, the Redfield tensor components by using coherent 
states can exphcitly be expressed as 

^xu^^n = 2 ^^1 ^^ 1^) ^^1 ^^ I'*) Ji^K^^){l + n{uj^f,)), if u^^ > 0, 

^xuf^n = 2 ^^1 ^z 1^) ^^1 ^^ I'*) J{^K^^)n{uJ^^), if u^^ > 0, (15) 

^X,y,^K = 2 ^^1 ^^ 1^^ ^'"l ^^ 1'*^ J(c^Ai;)(l + n{uJxv)), if UJXv > 0, 

Tai/^k = 2 ^-^l*^^ 1^^ ^'"I'^^l'*'^ '^(^''^)^*^^''^)' ifc^i,A>0, (16) 

where cu^ = Ef^/H are the system eigenfrequencies, n{uj) = i/(e^/''-^ — 1) is the distribution 
function for bosons, k, is the Boltzmann constant, T is the temperature, and J{uj) is the 
bath spectral function. 

To measure effects of the decoherence, one can use the entropy, the first entropy, and 
many other measures, such as maximal deviation norm, etc.( for example, see Refs. 23, 24, 
25). However, the decoherence of an open quantum system is essentially reflected through 
the decays of the off-diagonal coherent terms of its RDM. In general, the decoherence is 
produced due to the interaction of the quantum system with other system which has a 
large number of degrees of freedom, such as the devices of the measurement or environment. 
Here, we investigate the decoherence times via directly describing the evolutions of the off- 
diagonal coherent terms instead of using any measure of decoherence. In the following, we 
set the initial state of the qubit to p(0) = |(|0) + |1))((0| + (1|), which is a pure state and 
it has the maximum coherent terms, and the initial state of the environment is PhaMo) — 
Uke-f^^^/Tr{e-^^^^), where M^ = ookblbk, (3 = 1/kT. According to Ref. [12], we set uja = 
0.02 (ps)~^, Tc = O.luji in the calculations. 

We then shall use the Redfield equation to investigate the decoherence time of the DQD 
charge qubit. The analytical expressions of the elements of the RDM we obtained are 

^''^' l + 2n(cu2i) 2(l + 2n(cu2i))' ^^ 

^''^^^ ~ l + 2n(c.2i) + 2(l + 2n(c.2i))' ^^' 



_ (x + Vx' - culi)e(->^+V^'^^)* - (x - Vx' - cuii)e(->^-V^'^^)* 



+^UJ21 , , (19) 



P2l{t) 



WX^ - ^21 

-^^21 , ^^ ^ , (20) 

WX^ - ^21 

where x = ^-'^(^21) (1 + 2n(c<j2i))/2. Note that by substituting the parameters 0721, ujd, uji, 
ujc, 1] and T, in PCPB, DCPB and Ohmic bath into the above expressions, we can readily 
plot the evolutions of the p^2 (o^ P2i)- 

The evolutions of the coherent elements of the RDM of the DQD charge qubit in PCPB 
and DCPB with different Ui are plotted in Figs. 1 and 2 respectively. And when the 
environment is modeled with Ohmic bath the evolution of the coherent elements of the 
RDM is plotted in Fig. 3 as damping strength 77 = 0.04, 0.08 and 0.12. 

Ftg.l,Ftg.2,Ftg.3. 

It is shown that the decoherence times of the qubit decreases as the ui increases in PCPB and 
DCPB models. And the decoherence time of the qubit decreases as the damping strength 
f] in the Ohmic bath model increases. From Figs. 1-3, we see that the decoherence times 
of the DQD charge qubit in PCPB and DCPB are much shorter than the experimentally 
suggested result, which is in agreement with the results of Refs.[13, 15]. However, if we 
choose the environment as the Ohmic bath, the decoherence time of the qubit agrees with 
the experimental result when we choose a suitable coupling coefficient of the qubit to the 
environment. Moreover, we also consider the evolutions of the coherent term at different 
temperatures when the bath is modeled by PCPB, DCPB and Ohmic bath separately, which 
are plotted in Figs. 4, 5. and 6, respectively. 

Fig A, Fig.5, Fig.6. 

It is shown that the decoherence times decreases as the temperature increases no matter 
which bath is used to model the environment when ui (in PCPB and DCPB) and rj (in 
Ohmic bath) are fixed, this is physically reasonable. 



IV. CONCLUSIONS 

In summary, we have investigated the decoherence times of a double quantum dot (DQD) 
charge qubit when it is coupled to different baths by using the Redfield equation method. 
Our results show that the qubits have shorter decoherence times than the experimental ones 
as the environment is modeled by the acoustic phonon baths, which agrees with previous 
reports. Moreover, when we use Ohmic bath model the environment with the coupling 
coefficient of the qubit to the environment properly chosen, the decoherence time of the 
qubit is in well agreement with the experimental result. 
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V. FIGURES CAPTIONS 

Fig. 1: Pi2{t) vs time t. The evolution of the off-diagonal elements of the RDM for the 
DQD charge qubit in PCPB when ui = 0.5 (ps)~^ and uji = 0.7 (ps)^^. Here, Ud = 0.02 
(ps)^^, Qpz = 0.035 (ps)^^, and T = 30 mK. The initial state is described in the text. 

Fig. 2: Pi2{t) vs time t. The evolution of the off-diagonal elements of the RDM for the 
DQD charge qubit in DCPB. Here, Ud = 0.02 (ps)^^, Qdj = 0.029 (ps)^^. Other parameters 
are the same as those in Fig. 1. 

Fig. 3: Pi2{t) vs time t. The evolution of the off-diagonal elements of the RDM for the 
DQD charge qubit in Ohmic case, with cut-off frequency Uc = 0.05 (ps)^^, T = 30 mK and 
1] = 0.04, 0.08, 0.12, respectively. The initial state is described in the text. 

Fig. 4: Pi2{t) vs time t. The evolution of the off-diagonal elements of the RDM for the 
DQD charge qubit in PCPB when uji = 0.5 (ps)-\ T = 30 mK, 200 mK, 300 mK, 1 K 
respectively. Other parameters are the same as those in Fig. 1. 

Fig. 5: Pi2{t) vs time t. The evolution of the off-diagonal elements of RDM for the DQD 
charge qubit in DCPB when uoi = 0.5 (ps)~^, T = 30 mK, 200 mK, 300 mK, 1 K respectively. 



Other parameters are the same as those in Fig. 2. 

Fig. 6: Pi2{t) vs time t. The evolution of the off-diagonal elements of the RDM for the 
DQD charge qubit in Ohmic case when rj = 0.04 and T = 20 mK, 30 niK, 40 mK, 50 mK, 
respectively. Other parameters are the same as those in Fig. 3. 
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